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Abstract 

The U{1)b-l extension of the Standard Model requires the ex¬ 
istence of right-handed neutrinos and naturally realizes the seesaw 
mechanism of neutrino mass generation. We study the possibility of 
explaining the dark matter in this model with an additional scalar 
field, 4>dm, that is a singlet of the Standard Model but charged under 
U{1)b-l- An advantage of this scenario is that the stability of (poM 
can be guaranteed by appropriately choosing its B — L charge, with¬ 
out the need of an extra ad hoc discrete symmetry. We investigate in 
detail the dark matter phenomenology of this model. We show that 
the observed dark matter density can be obtained via gauge or scalar 
interactions, and that semi-annihilations could play an important role 
in the latter case. The regions consistent with the dark matter den¬ 
sity are determined in each instance and the prospects for detection 
in future experiments are analyzed. If dark matter annihilations are 
controlled by the B — L gauge interaction, the mass of the dark matter 
particle should lie below 5 TeV and its direct detection cross section 
can be easily probed by XENONIT; if instead they are controlled by 
scalar interactions, the dark matter mass can be much larger and the 
detection prospects are less certain. Finally, we show that this sce¬ 
nario can be readily extended to accommodate multiple dark matter 
particles. 
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1 Introduction 


About 27% of the energy-density of the Universe consists of an exotic form 
of matter commonly known as dark matter [1]. Throughout the years, many 
different dark matter candidates have been considered in the literature but 
so far no evidence has been found for any single one. Among them, Weakly 
Interacting Massive Particles (WIMPs) are particularly attractive for several 
reasons. First, particles with masses between the GeV and the TeV scale and 
with weak-strength interactions appear in many extensions of the Standard 
Model including the MSSM and UED models. Second, such particles natu¬ 
rally achieve, via a freeze-out in the early Universe, a relic density comparable 
to the observed value of the dark matter density -the so-called WIMP mira¬ 
cle. Third, WIMPs can be probed with current experiments in several ways, 
including direct [2] and indirect [3] dark matter searches as well as collider 
searches at the LHC |1]. In fact, we seem to be living through a very special 
time, the epoch in which the WIMP paradigm, which has dominated dark 
matter research for many years, is either going to be excluded or conhrmed 
0 . 

A problem common to many WIMP models is the explanation of the dark 
matter stability. It is indeed strange to have such a heavy particle (M ~ TeV) 
to be stable. Generally, dark matter stability is guaranteed by some discrete 
symmetry introduced exclusively for that purpose, usually a Z 2 . But that use 
of discrete symmetries is questionable not only due to its lack of motivation 
but also because they are expected to be broken by gravitational effects at the 
Planck scale, inducing dark matter decay and likely destroying the feasibility 
of such models. That is why it is often implicitly assumed that such discrete 
symmetries are actually the remnants of additional gauge |6] or flavor 
symmetries present at a higher scale, thereby delegating the problem to a 
framework larger than the model under consideration. It is not always clear, 
however, whether those assumptions are actually viable in such models. A 
better approach may be to hnd complete models where the dark matter is 
automatically stable due to the model’s structure and the quantum numbers 
of the dark matter held, without additional discrete symmetries W- Several 
models of this type have been studied in the literature [iniiiiiin], and they 
generally feature extended gauge sectors. 

In this paper, we consider a minimal extension of the Standard Model 
by a U{1)b-l gauge group. This extension is particularly appealing because 
it is remarkably simple, it is anomaly free once a right-handed neutrino per 
generation is introduced, and it naturally realizes the seesaw mechanism of 
neutrino mass generation. In addition, the U{1)b-l model has a rich collider 
phenomenology na dne to the presence of another neutral gauge boson. 
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Zbli that couples to quarks and leptons. This minimal setup, however, 
cannot account for the dark matter. To do so one needs either to modify the 
symmetry or increase the particle content. The former approach was applied 
in several works [HIISIIIHIIIT], where an additional symmetry was used 
to distinguish one of the particles, the would-be dark matter candidate. More 
natural seems the idea of increasing the particle content with an additional 
fermion or scalar field charged under B — L, with this charge chosen in such a 
way to guarantee its stability. For fermions, a vector-like pair is needed to not 
spoil anomaly cancellation, and the resulting scenario is very predictive, as 
emphasized recently in [IH]. Scalars, on the other hand, are not constrained 
by anomaly cancellation, so one more held is enough and, having both scalar 
and gauge interactions, typically give rise to a richer phenomenology. To our 
knowledge, scalar dark matter in this setup has been previously studied only 
in [12], which exclusively considered the scale invariant version of the B — L 
model. Here we want to study instead scalar dark matter within the general 
B — L extension of the Standard Model. Specihcally, we will analyze the 
different ways in which the relic density constraint can be satished, including 
resonant and non-resonant annihilations as well as semi-annihilations; we will 
determine the viable parameter space in each case; and we will investigate the 
detection prospects in future experiments. As we will show, when the dark 
matter density is determined by gauge-mediated interactions, the resulting 
framework is quite predictive and can be entirely probed by XENONIT. 
Finally, we also point out that multi-component dark matter can easily be 
accommodated within this scenario. 

The rest of the paper is organized as follows. In the next section we intro¬ 
duce our notation and describe in some detail the model. Our main results 
are presented in Section In it we numerically study the dark matter phe¬ 
nomenology of the model according to the interactions (gauge or scalar) and 
processes (annihilations or semi-annihilations) that set the relic density, and 
we analyze the detection prospects in each case. In Sectionj^we qualitatively 
discuss how to extrapolate our results to scenarios with multiple dark matter 
particles. Finally, we summarize our findings and draw our conclusions in 
Section [H 

2 The model 

Under Lf{l)B-L the Standard Model quarks and leptons have charge 1/3 
and —1 respectively. As a result, H — L is an anomalous symmetry and 
it is necessary to add additional fermions to consistently gauge it. In this 
paper, we consider a model based on the gauge symmetry SLf (3) x SLf (2) x 
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U (l)y X f/(l) b-l and containing three generations of right-handed nentrinos 
(iVj, i = 1,2,3) to cancel the gange anomalies, and a scalar, Sbli singlet of 
the Standard Model bnt charged nnder U{1)b-l-i that spontaneously breaks 
the B — L symmetry. We will assign charge -|-2 to S'bl so that the W acquire 
Majorana masses upon the breaking of the B — L symmetry, giving rise to a 
seesaw mechanism of neutrino mass generation. The Lagrangian then reads 

^ = ^SM- V{H, Sbl) - (^^XN^BBL^Ni + + h.c)j , (1) 

where the kinetic terms of the new particles have been omitted and, without 
loss of generality, it is assumed that the W mass matrix is diagonal after the 
spontaneous breaking of B — L. In the above equation, H is the Standard 
Model Higgs doublet and V {H, Sbl) includes all terms in the scalar potential 
involving Sbl (and possibly H too): 

V(H, Sbl) = i4S'blSbl - (S'blSbl) + y (2) 

In the unitary gauge H and Sbl can be written as 


H = 


Sbl — 


0 

h+VEW 

x/2 

hsL + VbL 

V2 - 


( 3 ) 

( 4 ) 


where vew = 246 GeV while vbl sets the B — L breaking scale, which can be 
much higher than the electroweak scale. The scalars h and Hbl inix with each 
other via the Xus term in equation (|^, giving rise to two mass eigenstates. 
Hi and H 2 , dehned by 


Hi = h cos 6 + Hbl (5) 

H 2 =—h sinO + Kbl cosO, (6) 

where 9 is the mixing angle. To satisfy the LHC bounds on the properties of 
the Higgs boson [201, 6 should be small. Hence, to a good approximation we 
can identify Hi with h and H 2 with Hbl- The scalar observed at the LHC 
with a mass of about 125 GeV [ 21 ] is then Hi. The mass of H 2 and 9 are 
instead free parameters of the model. 

The B — L gauge boson acquires a mass, after the spontaneous 

breaking of the B — L symmetry by Sbl, 

X^Zbl = ‘^QblVbl, (7) 
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being gsL the B — L gauge coupling constant. For simplicity, we assume a 
negligible kinetic mixing between U{1)b-l and f/(l)y so that Zbl and the 
SM do not mix with each other. The mass of Zbl and the gauge coupling 
qbl can be constrained with collider data. From LEP II data the bound 

My 

> 6 - 7 TeV ( 8 ) 

Qbl 

was derived in [22l|23]. At the LHC, current searches for dilepton resonances 
by ATLAS and CMS can also be used to set a bound on Zbl through the 
Drell-Yan process {qq Zbl ^ with i = e, fi) |21]. This bound turns 
out to be sligthly more stringent than the LEP II bound for Mzgj^ > 3 TeV 
-see [19]. 

After the breaking of the B — L symmetry, the right-handed neutrinos also 
become massive, with M^r- = Xj^-Vbl /(For simplicity, in our numerical 
analysis we take the right-handed neutrinos to be degenerate but our results 
hardly depend on this assumption). Thus, this model can naturally account 
for the observed pattern of neutrino masses and mixing angles via the seesaw 
mechanism. 

This minimal U{1)b-l extension, however, does not contain any dark 
matter candidates. We will therefore add one additional complex scalar held, 
4’ DM I to play that role. 4dm is a singlet of the Standard Model, has charge 
noM under U{1)b-l and does not acquire a vacuum expectation value. As 
we show below, tidm can be easily chosen in such a way that it is not possible 
to write a gauge invariant Lagrangian term that allows 4dm to decay, ren¬ 
dering 4dm automatically stable. In other words, the B — L symmetry itself 
stabilizes the dark matter particle, making unnecessary an ad hoc discrete 
symmetry. Regarding the dark matter phenomenology, 4dm will behave as 
a standard WIMP, obtaining its relic density via a freeze-out process in the 
early Universe, and giving rise to the usual signals in direct or indirect de¬ 
tection experiments, as we will describe in the next section. 

4dm has both gauge {B — L) and scalar interactions. The gauge interac¬ 
tions are determined, besides qbl and Mzg^, by udm and The scalar 

potential, on the other hand, now includes additional terms involving the 
dark matter held: 

V{4dm, H, Sbl) = Ldm4^dm4dm + Xh{4^dm4dm){H'^H) 

+ Xbm {4^Dm4dm){SbL^Bl) + -^4 (4 DM ^m) . (9) 

with Xh, Xdm and A 4 free parameters. The dark matter mass, 
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is then given by 


^Idm - ^ 


H- 


^EW 


+ A 


DM- 


■'BL 


> 0 . 


( 10 ) 


Notice that the dark matter mass, which is expected to be around the TeV 
scale to explain the observed dark matter density, receives a contribution 
proportional to vbl- If vbl were much larger than the TeV scale, say 
10^ GeV, an unnatural cancellation, or hne-tuning, between the hrst and 


third terms in equation (10) would be required to obtain the correct dark 


matter mass. To avoid such excessive hne-tuning, in our numerical analysis 
we take the dark matter mass as a free parameter but subject to the con¬ 
straint > '/XdmVbl/^- Even if somewhat arbitrary, this condition 

allows us to easily exclude the hne-tuned regions while scanning the param¬ 
eter space of the model. Thus, once we include scalar dark matter in the 
B — L model, the breaking scale, vbL) cannot be arbitrarily large and should 
lie relatively close to the TeV scale. In consequence, this setup gives rise to 
a low scale realization of the seesaw mechanism. 

Let us now discuss precisely what values can udm take to ensure dark 
matter stability within this model. Terms that would allow dark matter to 
decay are either Yukawa interactions of the form (j^DNiff or the scalar inter¬ 
actions (PdmHi, (pDMHiHj and (()DMHiHjHk, where Hi denotes either the SM 
Higgs doublet or Sbl- To prevent dark matter decay we must therefore choose 
riDM in such a way that all these terms are forbidden. Regarding the Yukawa 
interactions, (pDAiff, since (poM is a singlet under the SM gauge group, the 
combination //' should also be a singlet of SU{3)c x SU{2)b x U{1)y- And 
the only combination of fermions that produces such a singlet is / = N^, 
/' = Ni, yielding a term with B — L = ±2. Thus, for tidm ^ ±2 Yukawa 
interactions do not induce dark matter decay. Regarding scalar interactions, 
the mixing term (i )dm Sbl can be excluded with hdm 7^ ±2. Forbidding the 
cubic term (1)dmH^H requires udm 7^ 0, which means that cannot be 
a singlet oi B — L. From the other possible cubic term, ((>dmSbi, we get 
instead the condition udm 7^ ±4. Similarly, the quartic interaction (1)dmS'bb 
can be prevented with udm 7^ ±6. Finally, one may even forbid decays 
via dimension-5 operators such as 4>dmS%j^ (presumably suppressed by the 
Planck scale), which can be achieved for udm 7^ i8. Summarizing, the dark 
matter particle in this model will be stable provided that hdm 7^ ±2n, for 
n E Z and n < 4. Therefore, most choices of udm actually lead to a stable 
dark matter particle. 
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Figure 1: Left: The Feynman diagram corresponding to dark matter annihilation 
into fermions, fxjf —>■ ff, via the B-L gauge interaction. Right: The Feynman 
diagram for the spin-independent interaction between the dark matter and quarks. 


3 Scalar Dark Matter 

In this section we study the dark matter phenomenology of this model, in 
particular the (foM relic density and its detection prospects. Specihcally, we 
determine the regions in the parameter space of the model where the dark 
matter constraint can be satished according to the way in which the relic 
density is achieved in the early Universe. First, we will analyze the B — L 
gauge interactions and later the scalar interactions. 

3.1 Gauge interactions 

Since 0 dm cannot be a singlet oi B — L, it must necessarily have U{1)b-l 
gauge interactions. Let us then hrst consider the effects of these interactions 
on the dark matter phenomenology. Surprisingly, these effects have not been 
considered before, likely as a result of the strong bounds on qel and Mzg^. 
But as we will show, the B — L gauge interactions may not only explain 
the relic density but give also rise to a very predictive scenario; XENONIT 
should probe this possibility in the near future. 

The gauge interactions allow the annihilation of the dark matter particle 
into fermions mediated by the B — L gauge boson, (pff —)■ —)■ ff (see 

hgure[^ left), as well as the direct annihilation into two gauge bosons, 00'*^ —>■ 
ZblZbl- The relevant parameters in this case are then just four: the dark 
matter mass the B — L quantum number of 0dm (noM), the B — L 

gauge coupling (gBL), and the mass of the B — L gauge boson As 

we will see, due to the strong experimental constraints on Mzg^fQ bl, only 
the annihilation into fermions turns out to be relevant, and only close to the 
Zbl resonance. 

Dark matter annihilations are determined by the diagram shown in the 
left panel of hgure 1. Due to the structure of the scalar gauge coupling, 
the resnlting (av) is velocity suppressed (oc u^). All our numerical results 
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Figure 2: The dark matter relic density as a function of the dark matter mass for 
noM = 1/2 and three different values of (Mzg^, Qbl)- In this figure it is assumed 
that the gauge interactions determine the relic density. The horizontal (cyan) band 
displays the region consistent with cosmological observations. 


were obtained by implementing this model into micrOMEGAs [25] via Lan- 
HEP [2S]- Figure]^ shows the relic density as a function of the dark matter 
mass for n^M = 1/2 and three different values of {Mzgg, Qel), all of them 
compatible with collider bounds. The behaviour is analogous in all three 
cases, with the relic density decreasing as gets close to the resonance 

rs-/ Mzbl/‘^)j where it reaches its minimum value, and increasing as 
the dark matter mass moves away from it. Notice that the relic density tends 
to be much higher than the observed value (the horizontal cyan band) except 
in a narrow region close to the resonance. And even at the resonance, the 
relic density can be too large to be in agreement with the data, as illustrated 
by the solid (green) line. That is, sitting at the resonance may not guarantee 
a sufficiently small relic density. For the dash-dotted line, the impact of the 
annihilation into ZblZbl can be observed, as it gives rise to a slight sup¬ 
pression of the relic density once this channel becomes kinematically available 
2 TeV), but it is too small to have a signihcant effect. That is why 
only the annihilation into fermions mediated by Zbl is relevant in this case. 
From the hgure it is clear that, as expected, the higher the dark matter mass 
(or the gauge boson mass) the more difficult it is to satisfy the relic density 
constraint, and that for tidm = 1/2 it is not possible to satisfy it for dark 
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Figure 3: The viable parameter space when the dark matter relic density is deter¬ 
mined by the B — L gauge interaction. Left: Viable points projected onto the plane 
Mzgi^) for different values ofuDM^ 1/3 (green), 1/2 (red), and 1 (blue). 
Right: Viable points projected onto the plane gsb) for different values of 

udm- 


matter masses above 2 TeV or so. Next, we will investigate exactly what this 
upper limit on is and how it depends on tidm- 

To that end, we have randomly varied qbl and Mzgg for udm = 

1/2, 1/3, 1, and obtained a sample of viable models -those satisfying the relic 
density constraint, collider bounds and the perturbativity condition qbl < 
1. In the following we will analyze such a sample. Figure projects this 
set of viable models onto the planes Mzgg) in the left panel and 

qbl) in the right panel. The color convention is green, red and blue 
respectively for tidm = 1/3,1/2,1. From the left panel, we see that all the 
viable models indeed lie close to the resonance, with some non-negligible 
spread for n^M = 1- The range of dark matter masses increases with udm, 
extending from about 1 TeV for ubm = 1/3, to 2 TeV for hbm = 1/2, 
and to 4.5 TeV for ubm = 1- From the right panel we see that larger 
values of qbl can be obtained as udm increases. They reach 0.3, 0.6 and 
1.0 (the perturbative limit we imposed) respectively for ubm = 1/3,1/2,1. 
This hgure also explains why we hnd no viable models featuring > 

4.5 TeV: they require a gauge coupling larger than one -a possibility we did 
not consider. 

The same B — L gauge interaction that determines the relic density in¬ 
duces a coupling between the dark matter and quarks that can be probed 
in direct detection experiments -see hgure right. Being the dark matter 
candidate a scalar, only the spin-independent cross section, asi, is relevant. 
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Figure 4: The spin-independent direct detection cross section as a function of the 
dark matter mass for our set of viable points, classified according to the value of 
noM- 1/3 (green), 1/2 (red), and 1 (blue). All these points satisfy the relic density 
constraint via gauge-mediated dark matter annihilations. For comparison, we also 
display the current bound from LUX (solid line) and the expected sensitivity of 
XENON IT (dashed line). 


Figure 1^ displays our sample of viable models in the plane a si) us¬ 

ing the same color convention as before. For comparison, the present bounds 
from XENONIOO [27] and LUX [2S] (solid lines) are also shown, as well as 
the expected sensitivity of XENONIT [29] (dashed line). When udm = 1/3 
(green points), the models lie in a tiny area featuring < 1 TeV and 

2.7 < (T5//10“®pb < 5 -very close to the LUX bound. A small fraction of 
them are in fact already excluded and the rest can be easily probed in future 
experiments. When uom = 1/2 (red points), the models lie in a larger area 
extending to dark matter masses of order 2 TeV and with cross sections be¬ 
tween 10“® and 10“® pb. This region will be entirely probed by XENONIT. 
Finally, when tidm = 1 (blue points), the models occupy a broader area 
that extends down to cross sections of order 3 x 10“^° pb. Practically all 
of them lie within the expected sensitivity of XENONIT. And the few lying 
below the XENONIT sensitivy can be probed by either XENON-nT or LZ 
[30] . Thus, this scenario can be tested in the near future via direct detection 
experiments. 

Regarding indirect detection, no signals are expected in this case because 
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the annihilation rate today is velocity snppressed and therefore negligibly 
small. This fact conld be nsed to falsify this scenario if a credible indirect 
detection signal of dark matter were observed. 

Snmmarizing, dark matter annihilations mediated by the B — L gange 
boson allow to satisfy the relic density constrained bnt only close to the res¬ 
onance ~ Mzgj^) and within a narrow range of dark matter masses 

that depends on udm- Moreover the spin-independent direct detection cross 
section is expected to be within the reach of the XENONIT experiment, 
providing a direct way to probe this scenario. 

3.2 Scalar interactions 

Besides the gange interactions we have already stndied, the dark matter 
particle, (poM, also has scalar interactions with both H ^ Hi (the Higgs) and 
Sbl ~ H 2 -see eqnation (|^. Their effect on the dark matter phenomenology 
was partially stndied in [T9|, bnt only within a scale invariant framework. 
Here we generalize their findings. 

Scalar interactions between the dark matter and the SM scalar give rise 
to the well-known Higgs-portal scenario [211 ESI E3|, whose phenomenology 
has been thoronghly stndied in the recent literatnre -see e.g. [MlEHlEniEIl 
ESI ESj. Notice however that when the scalar Higgs-portal arises from the 
B — L model, one antomatically obtains an explanation for the stability of the 
dark matter not based on an ad hoc Z 2 symmetry. This fact may be relevant 
in view of the possible effects dne to dark matter decay indnced by the 
breaking of discrete symmetries at the Planck scale, as recently emphasized 
in BDl- In addition, the dark matter field is necessarily complex -charged 
nnder U{1)b-l- rather than real, as nsnally assnmed in standard Higgs- 
portal models. 

In the following we will focns on the interactions with H 2 or, to a good 
approximation, with Sbli the scalar that breaks the B — L symmetry. Snch 
interactions give rise to two distinct hnal states in dark matter annihilations: 
NiNi (mediated by H 2 ) and H 2 H 2 . The relevant parameters in this case are 
the right-handed nentrino masses (MjvJ, the mass of H 2 the scalar 

conpling {Xdm), the vev of Sbl (vbl), and the dark matter mass 

Fignre shows the relic density as a fnnction of the dark matter mass 
for \dm = 0.2 and different valnes of Mh 2 '- 1 TeV (dash-dotted line), 2 TeV 
(dashed line), and 3 TeV (solid line). The common right-handed nentrino 
mass was taken to be 0.4 TeV. The effect of the H 2 resonance is clearly 
observed in this figure as it leads to a strong suppression of the relic density 
for ~ Mji/ 2 / 2 , reaching values as low as 10“^. In fact, for the chosen 

value of Xdm, the relic density constraint (cyan horizontal band) can only 
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Figure 5: The dark matter relic density as a function of the dark matter mass 
for Xdm = 0.2 and different values of Mh 2 - 1 TeV (dash-dotted line), 2 TeV 
(dashed line), and 3 TeV (solid line). The horizontal (eyan) band displays the 
region consistent with cosmological observations. 


be satisfied close to the H 2 resonance. The effect of the other annihilation 
channel, H 2 H 2 , can also be observed for Mh 2 = 1 TeV (dash-dotted line), as 
it induces a slight decrease in the relic density when > 1 TeV. 

To illustrate the dependence with Xdm, we show in figure the relic 
density as a function of for Mh2 = 1 TeV and Xdm = 0.2 (dash- 

dotted line), 0.4 (dashed line), and 0.6 (solid line). As expected, the larger 
the coupling, the smaller the relic density. A qualitative difference is that 
the larger couplings allow to satisfy the relic density constraint far from the 
H 2 resonance, via the annihilation into H 2 H 2 . For Xdm = 0.4, the relic 
density is below the observed value for dark matter masses between 1 TeV 
and 1.3 TeV while for Xdm = 0.6 this happens for between 0.9 TeV 

and 1.7 TeV. In the following we will study in more detail the resulting 
viable parameter space. To do so, it is convenient to distinguish between the 
two possible final states. First, we will address dark matter annihilation into 
right-handed neutrinos {(fcff —)■ NiNf), which, as we have seen, is relevant 
close to the H 2 resonance (2M<^^^ ~ Mh 2 )- Second, we discuss dark matter 
annihilation into H 2 H 2 , which is possible for Mh 2 < ^iPdm- 

Let us start with <Pdm<P^dm which proceeds through an s-channel 

diagram mediated by H 2 , giving rise to a (av) proportional to {XdmVblY■ 
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Figure 6: The dark matter relic density as a function of the dark matter mass 
for Mh2 = 1 TeV and different values of Xdm- 0.2 (dash-dotted line), 0.4 (dashed 
line), andO.6 (solid line). The horizontal (cyan) band displays the region consistent 
with cosmological observations. 



Figure 7: The viable parameter space when the relic density is determined by dark 
matter annihilation into right-handed neutrinos XhiXfi- Left: The 

viable points projected onto the plane (M^p^^, Mh^). Right: The viable points 
projected onto the plane Xdmvbl)- 


Interestingly, this annihilation process involves three of the four new helds 
included in this U{1)b-l extension -only Zbl is missing. In [12] this process 


13 









4’\ z' tl2 

\ / 

^ H 2 ^ 

> - ( 

/ \ 

/ \ 


■H. 


■ Ho 


Hi 

\ X ^ 

N X 

\ y 

X 

/ \ 

/ N 

V S.H2 


Figure 8: The Feynman diagrams that contribute to the annihilation of dark mat¬ 
ter into the final state H 2 H 2 . 


was briefly mentioned but not analyzed. As before, we have scanned the rel¬ 
evant parameter space and selected the points compatible with accelerator 
and perturbativity bounds, and with the dark matter constraint, requiring 
the relic density to be determined by the annihilation into right-handed neu¬ 
trinos. The resulting viable models are shown in £gure[^ projected onto the 
planes Mh^) in the left panel, and XdmHbl) in fhe right 

panel. From the left panel we see that for < 2 TeV the viable models 

lie along two branches, one slightly above and the other slightly below the H 2 
resonance. It is difficult to hnd viable models right on top of the resonance 
because the relic density tends to be quite suppressed there, as we already 
saw in hgures and But as the dark matter increases, the two branches 
tend to merge and the viable models move closer and closer toward the reso¬ 
nance. At very high masses, the relic density constrained is satished only on 
top of the resonance. Regarding the dark matter mass range, we see that, 
as a result of the resonant annihilation, it can extend to values as large as 
50 TeV. From the left panel we see that, as expected from the cross section, 
the product XdmVbl must increase with the dark matter mass in order to 
maintain the correct value of the relic density, reaching, for ~ 50 TeV, 

a value of order 100 TeV. For vbl we hnd in our sample a maximum value 
of about 200 TeV. We see then how the dark matter constraint can set an 
upper bound on the otherwise unconstrained B — L breaking scale. 

Now let us move to the second case, H 2 H 2 . This is a non¬ 
resonant process requiring Mh 2 < that can proceed through an s- 

channel diagram mediated by i/ 2 , a t- and u-channels mediated by 0 _dm, 
and via direct annihilation -see hgure This process was considered in 
detail in [19], where some analytical formulas can be found. Their numerical 
results, however, strongly depend on the assumed scale invariance and differ 
significantly from ours, as illustrated below. 

We display the viable models in the plane Xdm) in the left panel 

of hgure As expected, the minimum value of Xdm increases with the dark 
matter mass so as to maintain the correct relic density. The maximum value 
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Dark Matter Mass [TeV] Dark Matter Mass [TeV] 

Figure 9: The viable parameter space when the relic density is determined by dark 
matter annihilation into the heavier scalars -(t>DM4>^£iM H 2 H 2 . Left: The viable 

points projected onto the plane Xdm)- Right: The viable points projected 

onto the plane 


of the dark matter mass in this case lies around 9 TeV, where the observed 
dark matter density can be obtained only for Xdm close to 1 -the perturbative 
bound we have imposed. In the scale invariant version of this model, studied 
in [in], it was instead found that the relic density constraint could only be 
satished for Xdm ^ 1- The right panel displays the viable models in the 
plane dark matter mass versus the product which controls 

the relative contribution between the direct annihilation amplitude and the t- 
and M-channels -see [m]. Since it is larger than 1, it is the t- and u-channels 
that dominate the annihilation rate, and they are more dominant at larger 
dark matter masses. In addition, that same product measures the hne-tuning 
required to obtain the dark matter mass, and we require it to be smaller than 
9 -see Section From the hgure we see that such an upper bound is often 
reached and that it helps set the upper limit on the dark matter mass. 

The interaction between the dark matter and quarks proceeds in this case 
via diagrams mediated by Hi, H 2 and Zbl, which are not constrained by the 
relic density -they depend on additional parameters such as 9, Xh, gsL, and 
Mzg ^. As a result, the direct detection cross section can vary over a very wide 
range. Regarding indirect detection, the advantage of the scalar interactions 
is that the annihilation rate is not velocity suppressed and it takes instead 
the expected value, an ~ 3 x 10“^®cm^s“^. But since current dark matter 
indirect detection experiments (mostly those searching for gamma rays jH] 
and antiprotons m ) only probe the small mass region (M < 100 GeV), no 
strong bounds can be obtained on this scenario. 
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Figure 10: The dark matter relic density as a function of the dark matter mass 
for Mh 2 = 200 GeV and different values of X3: 0.1 (dash-dotted line), 0.5 (dashed 
line), 1.0 (solid line). The horizontal (cyan) band displays the region consistent 
with cosmological observations. 


3.3 Semi-annihilations 


The scalar interactions described in the previous section are generic, they are 
present for any value of udm- One can also have, however, interactions which 
are allowed only for specific values of udm- For instance, when udm = 2/3 
the tern£] 


=2^2/3 — ~ 


Sbl + h.c. 


dm^bl 


( 11 ) 


becomes compatible with the B — L gauge symmetry. This term allows the 
dark matter particles to semi-annihilate |13] via the process (fcf —)■ (j)^H 2 
{(fxf) —)■ (f^Hi also contributes but is additionally suppressed by the mixing 
angle). For simplicity, we will assume in this section that all other interac¬ 
tions of the dark matter particle are irrelevant. Thus, the dark matter phe¬ 
nomenology is determined by only three parameters: A 3 , and Mh 2 - 

Far from the H 2 threshold we get the simple behavior (av) oc 

Figure 10 shows the relic density as a function of the dark matter mass 
for Mh 2 = 200 GeV and A 3 = 0.1, 0.5, 1.0. For comparison the region 
consistent with current observations is also displayed as a horizontal band. 


^Another possibility is udm = !> which renders the term gauge invariant 

but it does not gives rise to semi-annihilations. 
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Figure 11: The viable parameter space when the relic density is determined by 
semi-annihilations, projected onto the plane X^). 


The behavior of the relic density for a given value of A 3 is easy to understand: 
it simply increases with the square of the dark matter mass except close to 
the kinematic threshold ~ Mh 2 )- Notice that for A 3 = 0.1, the 

resulting relic density is always above the observed range. For A 3 = 0.5 and 
A 3 = 1 compatibility with the observed data can be obtained respectively 
for 700 GeV and ~ 1.3 XeV. Since we require A 3 < 1, we 

can already conclude that semi-annihilations do not allow to satisfy the relic 
density constraint for dark matter masses above 1.3 TeV or so. 

This fact is further illustrated in hgure IT, which shows a scatter plot 
of viable models in the plane A 3 ). For this hgure we assumed that 

semi-annihilations determine the relic density. Since A 3 should increase with 
the dark matter mass to keep the relic density constant, most points lie along 
a band in this plane. The upper bound on the dark matter mass is obtained 
when A 3 = 1 and it reaches about 1.4 TeV. 

Since the interaction in equation ( pT| cannot give rise to dark-matter 
nucleon scattering, there is no prediction for the spin-independent direct 
detection cross section, which will be determined by the usual diagrams me¬ 
diated by Hi 2 and Zbl- The annihilation rate today, on the other hand, 
does depend on A 3 . In fact, it was stated in [19] that this semi-annihilation 
process could account for the gamma ray excess from the Galactic center 
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claimed in jUj. 


4 Multi-component Dark Matter 

As we have seen in the previous sections, by simply adding an extra scalar 
held charged under B — L, the U{1)b-l extension of the Standard Model 
can account for the observed dark matter in the Universe. The stability of 
the dark matter is then guaranteed by the choice of the B — L dark matter 
charge whereas the observed relic density can be obtained via annihilations 
induced by gauge or scalar interactions, or via semi-annihilations. It turns 
out that this minimal framework can also be extended in a straightforward 
way to a dark matter sector consisting of multiple dark matter particles, a 
scenario known as multi-component dark matter. 

Many models, as a result of the discrete symmetry that is used to stabilize 
the dark matter, only allow for a single dark matter particle -typically the 
lightest odd state under a Z 2 symmetry. But there is currently no evidence 
that the dark matter sector is so simple. In principle, it could consists of 
a number of distinct particles with different properties, just like the visible 
sector. Models with multi-component dark matter have certainly been stud¬ 
ied before, particularly those with two-component dark matter. And they 
typically require larger discrete symmetries, such as Z 2 x Z 2 |l5lll6], to stabi¬ 
lize the dark matter particles (an exception is the model considered in |TF]h 
Within the B — L model, multi-component dark matter has the advantage of 
not requiring discrete symmetries at all and of allowing even for more than 
two dark matter particles. 

Let us qualitatively describe how this extension would work out for two 
dark matter particles, and (j) 2 - To avoid mixing between them, their B — L 
charges must be different, say 1/2 and 1/3. As before, this choice guaran¬ 
tees that both will automatically be stable without the need of an additional 
discrete symmetry. Regarding the relic density, there are essentially two dif¬ 
ferent ways in which it can be obtained. The hrst one is a combination of the 
possibilities mentioned in the previous section. For instance, 0i could anni¬ 
hilate resonantly via the gauge interaction whereas 02 annihilates via scalar 
interactions into H 2 H 2 ] or 0i could annihilate instead into right-handed neu¬ 
trinos at the H 2 resonance. The second way is via the annihilation of dark 
matter into dark matter; that is, the heavier dark matter particle annihilates 
into the lighter one (0202 —t 0i0i) via the new scalar interaction 



( 12 ) 
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which is always allowed. The lighter dark matter particle would then anni¬ 
hilate through one of the ways mentioned in the previous section. 

Even though this idea can in principle be extrapolated to any number of 
dark matter particles ( 01 ,..., 0 n), there is a limit on this number imposed by 
the dark matter constraint, which becomes stronger in models with multiple 
dark matter particles. As evidenced by the fact that it is often necessary to 
rely on resonances to satisfy the dark matter bound, the scalar relic density 
tends to be high in this model. If we now include several dark matter parti¬ 
cles, the observed dark matter density would be the sum of the relic densities 
of all of them. Thus, each particle should contribute less than required in 
the previous section, making it more difficult to satisfy the dark matter con¬ 
straint. But it is not hard to hnd viable regions for a small number of dark 
matter particles. 


5 Conclusions 

We studied in detail scalar dark matter in the B—L extension of the Standard 
Model. In it, the SM particle content is extended by three right-handed 
neutrinos required to cancel the gauge anomalies, one scalar held to break 
the B — L symmetry, and another scalar held charged under B — L to explain 
the dark matter. This model naturally incorporates the seesaw mechanism 
of neutrino mass generation while dark matter stability is easily achieved by 
appropiately choosing the B — L charge of Moreover, since the dark 

matter mass receives a contribution proportional to vbl, the relic density 
constraint implies that the scale oi B — L breaking cannot be much above 
the TeV scale. 

The dark matter particle in this scenario has both gauge and Yukawa 
interactions and its relic density can be the result of annihilations or semi¬ 
annihilations. We analyzed these diherent possibilities and determined the 
viable parameter space in each case. If the dark matter relic density is de¬ 
termined by the B — L gauge interactions ^bl —t //), the 

annihilations should happen very close to the resonance 
and the dark matter mass cannot be larger than about 4.5 TeV (for udm ^ 
1). In addition, the spin-independent cross section is expected to be large 
enough to be probed by the upcoming XENONIT experiment. Scalar inter¬ 
actions, on the other hand, induce annihilations into right-handed neutrinos 
( 0 DM 0 DM H 2 ^ NiNi) and into the heavier scalars {(j)DM4>DM —^ H 2 H 2 ) 
via different channels. The former process takes place resonantly and al¬ 
lows the dark matter to be as large as 50 TeV whereas the latter one is 
non-resonant and requires a dark matter mass below 9 TeV to be in agree- 
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ment with the observed value of the dark matter density. When udm = 2/3 
semi-annihilation processes t become possible and can 

account for relic density provided that ^1-4 TeV. Finally, we showed 

that one can easily accommodate multi-component dark matter in this sce¬ 
nario. The U{1)b-l extension of the Standard Model therefore provides a 
compelling and testable framework for neutrino masses and scalar dark mat¬ 
ter. 
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